Effect of pressure on the electronic and magnetic properties of CdV 2 04: Density 

functional theory studies 
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We investigate the effect of pressure on the electronic and magnetic states of CdV204 by using 
ab initio electronic structure calculations. The Coulomb correlation and spin-orbit coupling play 
important role in deciding the structural, electronic and magnetic properties of the compound. The 
total magnetic moment of V ion is found to be — 1.3 \ib and making an angle of —9.5 degree with 
the z-axis. In the tetragonal phase, the ground state is the orbital ordered state where V d xz 
and d yz obtitals are mainly occupied at the neighbouring sites. This work predicts the electronic 
phase transition from orbital-ordered-insulator to orbital-ordered-metal to orbital-disordered-metal 
with increasing pressure. The pressure induced broadening of lower and upper Hubbard bands 
gives rise to metal-insulator transition above 35 GPa. The simple mean-field theory used in the 
present work is able to describe the pressure dependent variation of the antiferromagnetic transition 
temperature suggesting the applicability of the method in the study of the magnetic behaviour of 
similar geometrically frustrated systems. 

PACS numbers: 75.25.Dk, 71.20.-b, 71.27.+a 



I. INTRODUCTION 

The spinel vanadates having general formula AV2O4 
(A=Cd, Mg and Zn) offer a fertile ground to study many 
aspects of condensed matter physics^— At room temper- 
ature they are paramagnetic insulator with face-centred 
cubic structure. The four V ions forming a regular tetra- 
hedron and their spins are antiferromagnetically coupled 
resulting in an interesting example of a geometrically 
frustrated system. At low temperature they show struc- 
tural transition from cubic to tetragonal and the tetra- 
hedron becomes distorted as shown in Fig. 1. The four 
sites of the V atoms forming the tetrahedron are or- 
bitally inequivalent as VI and V2 atoms have different 
orbital occupancies. The magnetic transition from para- 
magnetic to antiferromagnetic (AFM) is also observed at 
lower temperature 2^— 

Among the above mentioned three compounds the 
structural transition temperature (T$) is highest (—97 



K) for CdV 2 4 and lowest (-50 K) for ZnV 2 4 
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On 



the other hand the AFM transition temperature (T/v) 
is lowest (-35 K) for CdV 2 4 and highest (-42 K) 
for MgV 2 4 4 ' 12 i 14 The frustration index, defined as / 
= \9cw\/Tn where 8cw is the Curie- Weiss tempera- 
ture and TV is the AFM transition temperature, is con- 
sidered as a measure of the activeness of geometrical 
frustration^! 4 .^ The values of / for CdV 2 4 , MgV 2 4 , 
and ZnV 2 4 compounds are found to be about 11.4, 
14.3, and 21.3, respectively, suggesting that the geomet- 
rical frustration is least active in CdV 2 4 J^ The pres- 
sure dependent magnetization studies on CdV 2 4 and 
ZnV 2 4 have been carried out by Canosa et al. and 
found opposite trends for the T/v. The T/v of CdV 2 4 
and ZnV 2 4 , respectively, increases and decreases with 
increase in pressured These contrasting behaviours in 
this series of compounds are unusual and yet to be un- 



derstood. 

It is clear from the above discussion that CdV 2 4 is 
a relatively simple compound where the effect of spin 
fluctuations due to geometrical frustration is expected 
to be minimal. Thus CdV 2 4 would be a good start- 
ing candidate to study the complex physical properties 
of the above compounds using density functional theory 
(DFT), which is a mean- field theory. Here we employ 
the DFT based ab inito electronic structure calculations 
to study the pressure dependent electronic and magnetic 
properties of CdV 2 4 compound. The Coulomb correla- 
tion among V 3d electrons are found to be important 
in understanding its structural and electronic proper- 
ties. The spin-orbit coupling (SOC) is weak and total 
magnetic moment of the V ion comes out to be —1.3 
fiB which is closer to the experimental value of —1.2 
fiB- It shows orbital ordered ground state where d xz 
and dy Z orbitals are mainly occupied at the neighboring 
V sites. The electronic states of the compound evolves 
from orbital-ordered-insulator to orbital-ordered-metal 
to orbital-disordered-metal with increasing pressure. The 
broadening of lower and upper Hubbard bands due to 
increasing pressure is found to be responsible for the 
pressure induced metal-insulator transition. The present 
mean-field calculations correctly describe the pressure de- 
pendent shift of Tjv indicating the applicability of such 
approach in understanding the magnetic properties of 
similar systems. 



II. COMPUTATIONAL DETAILS 

The ferromagnetic (FM) and AFM calculations of 
CdV 2 4 have been carried out by using state-of-the- 
art full-potential linearized augmented plane wave (FP- 
LAPW) method. 17 In the AFM solution the directions of 
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spins of VI and V2 atoms shown in Fig. 1 are consid- 
ered to be opposite. The muffin-tin sphere radii used in 
the present work are 2.0, 1.85, and 1.52 Bohr for Cd, V, 
and O atoms, respectively. In order to see the effect of 
pressure on the electronic and magnetic properties of the 
compound only volume of the primitive unit cell is varied. 
The atomic positions and tetragonal distortion are kept 
fixed to the reported experimental values of 85 KJ£ For 
the exchange correlation functional, we used generalized 
gradient approximation (GGA) form of Perdew et alJ2- 
The effect of on-site Coulomb interaction among V 3d 
electrons is also considered within GGA+U formulation 
of the DFT 2& The values of U and J used in the calcula- 
tions are 4.0 and 0.5 eV, respectively, which were found to 
describe the correct electronic and magnetic properties of 
MgV2 0*4 compound^ The self-consistency was achieved 
by demanding the convergence of the total energy to be 
smaller than 10~ 5 Hartree/cell. To study the role of or- 
bital degree of freedom in deciding the magnetic prop- 
erties of the compound the spin-orbit coupling (SOC) is 
also considered. The initial axis of magnetization was set 
along the z-axis and convergent solution was obtained in 
presence of SOC, which corresponds to a non-collinear 
magnetism. The consideration of SOC in the calculation 
is very time consuming, therefore for this solution the 
energy convergence criterion was set to be smaller than 
10~ 4 Hartree/cell. 



III. RESULTS AND DISCUSSIONS 

The energy vs volume plots of CdV204 obtained from 
FM and AFM GGA+U calculations are show in Fig. 2. 
Both the curves show almost parabolic behaviour and 
the volume corresponds to minimum energy is found to 
be around 1092 Bohr 3 . The energy of AFM solution is 
always found to be less than that of FM solution suggest- 
ing the AFM ground state of the compound for the range 
of unit cell volumes studied here. The energy difference 
between FM and AFM phase increases with decrease in 
volume. Thus decrease in volume gives rise to more sta- 
bility to the AFM phase in comparison to FM phase. In 
order to know the exact equilibrium volume of the com- 
pound we fitted the energy-volume data by using equa- 
tion of states formula of Vinett et al^- The equilibrium 
volumes thus obtained correspond to FM and AFM so- 
lutions are 1096.6 and 1093.1 Bohr 3 , respectively, which 
are about 1% and 1.3% less than the experimental value 
of 1108.2 Bohr 3 at 85 The bulk moduli of the com- 

pound found for FM and AFM solutions are ~180.8 and 
~170.7 GPa which is about 0.4 times less than the bulk 
modulus of diamond. For U—5 eV the equilibrium vol- 
ume and bulk modulus of the ground state AFM phase 
is found to be about 1098.4 Bohr 3 and 175.5 GPa, re- 
spectively, suggesting that the increased value of U does 
not have much influence on its structural properties. The 
present calculated value of the bulk modulus is about 20 
GPa more than that reported earlier and obtained from 



the use of hybrid functional. 16 

The partial density of states (PDOS) of Cd 4d, V 3d 
and O 2p states obtained from FM GGA+U solution 
corresponding to equilibrium volume are shown in Fig. 
3. The valance band can be divided in three regions: (i) 
region A from to -2.0 eV, (ii) region B from -2.2 to -7.8 
eV and (iii) region C below -7.8 eV. Regions A and C 
are dominated by V 3d and Cd 4d states, respectively, 
whereas region B has similar contribution from all the 
three atoms. The shape of Cd Ad, V 3d and O 2p PDOS 
are similar in all the regions suggesting the hybridization 
between Cd 4d, V 3d and O 2p orbitals. The lowest 
lying conduction bands, denoted by 1 and 2 in Fig. 3(b), 
have dominating contribution from V 3d states. There 
is a finite gap between regions A and B which can be 
attributed to the gap between bonding and antibonding 
bands arising due to finite overlap of V 3d and O 2p 
orbitals. The band gap of ~0.3 eV is essentially arising 
due to correlation induced splitting of V 3d bands in lower 
and upper Hubbard bands putting this compound into 
the category of Mott-Hubbard insulator. 

In order to know the exact magnetic moment (MM) of 
the V ion we have included SOC in the calculation and 
obtained the FM solution. The spin part of MM at the 
V site is found to be about 1.6 \jlb and directed along the 
z-axis. The total spin only MM is found to be 2 \xb per 
formula unit (fu) suggesting that effectively V ion is in 
5=1 spin state as expected for 3+ ionic state of the V 
in this compound. The orbital moment is mainly lying 
in the yz-plane. The y and z components of the orbital 
moment (L) are found to be about -0.1 and -0.2 respec- 
tively where negative sign indicates the direction of the 
moment. Since the direction of the spin moment (S) is 
mainly along the z-axis. Thus the total moment (J) also 
lies in the yz-plane. The y and z components of J are 
found to be about -0.1 and 0.6. Using these values of L, 
S and J one can estimate the magnitude and direction of 
total magnetic moment of V ions. The magnitude comes 
out to be ~1.3 fJ-B, which is in good agreement with the 
experimental value of about 1.2 [is ^ The total magnetic 
moment is making an angle of ^9.5 degree with the z- 
axis. Here it is important to note that there is a large 
difference between the calculated and experimentally ob- 
served MM < ■: M . V ( ) . compound and this behaviour has 
been attributed to the activeness of geometrical frustra- 
tion responsible for quantum fluctuations However, 
the small difference of 0.1 /is between the calculated and 
experimental MM for CdV204 suggests that the geomet- 
rical frustration may be less important for understanding 
its magnetic properties. This observation can further be 
justified by comparing the frustration index / of both the 
compounds. The values of / for CdV204 and MgV204 
have been reported to be 11.4 and 14.3, respectively^ 
suggesting that geometrical frustration is less active in 
the CdV204. Finally one should keep in mind that the 
DFT is a mean-field theory and effect of quantum fluc- 
tuations can not be studied using this technique. There- 
fore, in the geometrically frustrated system the difference 



3 



between the experimental and calculated magnetic mo- 
ment is expected; and it can be directly related to the 
/. Since the value of / for MgV2C>4 is greater than that 
of CdV2C>4, therefore the difference between the experi- 
mental and calculated magnetic moment for MgV204 is 
expected to be more in comparison to CdV2C>4 as men- 
tioned above. 

The electronic occupancy of d xz (d yz ) orbital for two 
vanadium atoms VI and V2 forming the tetrahedron (see 
Fig. 1) is found to be about 0.79 and 0.06 (0.06 and 0.79), 
respectively. The occupancies of rest of three V 3d or- 
bitals are same at each V sites. This result clearly shows 
the orbital ordered ground state for the CdV2C>4 com- 
pound where d xz and d yz orbitals are mainly occupied at 
VI and V2 sites of the tetrahedron. Here it is important 
to note that the AFM solution does not have much in- 
fluence on the orbital ordering and solely arising due to 
tetragonal structure of the compound. 

The evolution of electronic structure of the compound 
with increasing pressure is shown in Fig. 4, where we 
have plotted total density of states (TDOS) per formula 
unit (fu) obtained from AFM solutions at different pres- 
sure. One can notice couple of obvious changes in the 
TDOS with varying pressure. With increasing pressure 
Cd 4d states are becoming deeper in energy. The gap 
between regions A and B enhances with increasing pres- 
sure suggesting the increased energy difference between 
the bonding and antibonding molecular orbitals mainly 
consist of V 3d and O 2p atomic orbitals. This behaviour 
is as per the expectation because increasing pressure in- 
creases the overlap between V 3d and O 2p atomic or- 
bitals and the energy difference between bonding and an- 
tibonding molecular orbitals is proportional to the over- 
lap between the atomic orbitals. With increase in pres- 
sure the band gap decreases from —0.5 to —0.27 eV at 
16.6 GPa as conduction band (CB) 1 comes closer to the 
valance band (VB). The gap almost vanishes at 34.5 GPa 
as the maximum of VB and minimum of CB 1 merges 
with each other as evident from Fig. 4(c). On further 
increase of pressure the VB and CB 1 started overlapping 
and such overlap of VB and CB may lead to metallic con- 
ductivity above 34.5 GPa. The TDOS/fu at the Fermi 
level increases from — 1 states/eV at 87.5 GPa to —2 at 
105.5 GPa enhancing the conductivity by 2 times. Dis- 
persion relations along high symmetric directions of the 
Brillouin zone for 0.23 and 34.5 GPa are shown in Fig. 
5. It is evident from Fig. 5(a) that the valance band 
maximum and conduction band minimum lie at the Z 
point and T point, respectively. The energy difference 
between these two points gives the indirect band gap 
of —0.5 eV, which is responsible for the insulating be- 
haviour of the compound. When pressure is increased 
to 34.5 GPa the energy difference between the Z and T 
points becomes almost zero. This suggests the pressure 
induced insulator to metal transition above 34.5 GPa. 
The increase in pressure brings V atoms closer resulting 
in increased overlap between V 3d orbitals. This finally 
leads to broadening of upper and lower Hubbard bands 



giving rise to reduced separation between them which is 
responsible for the decrement in the gap with increasing 
pressure. Thus present work clearly indicates the pres- 
sure induced metal-insulator transition in the compound 
due to broadening of upper and lower Hubbard bands. 

To study the effect of pressure on the orbital order- 
ing we have plotted the pressure dependent electronic 
occupancies of V d xz and d yz orbitals at VI sites ob- 
tained from AFM solutions in Fig. 6. The degree of 
orbital polarization at VI sites can be defined as P or j, 
= (dxz-dy Z ) / (d xz +d yz ). P orb = 1 corresponds to 100% 
orbital polarization as only d xz orbital is occupied at 
VI sites. Needless to say that this situation automat- 
ically tells about 100% orbital polarization at V2 sites 
where only d yz orbital is occupied. P or b = corre- 
sponds to orbital disordered state as same kind of or- 
bitals will be occupied at every sites. The occupancies 
of both the orbitals show almost linear pressure depen- 
dence in the insulating region and non-linear pressure 
dependence in the metallic region. The occupancy of d xz 
continuously decreases with increase in pressure. The 
pressure dependent occupancy of d yz orbital shows non- 
monotonous behaviour. It decreases in the insulating 
region and increases in the metallic region. The P orb 
also shows non-monotonous pressure dependence. In 
the insulating phase the value of P or b increases slowly 
from -0.46 at GPa to -0.52 at 34.5 GPa. However, 
the decrement of P or b is fast in the metallic phase as 
its value reduces from —0.51 at 71.6 GPa to —0.32 at 
105.5 GPa. This result suggests the complete melting 
of the orbital ordering at higher pressure. Thus, the 
present work clearly shows the interesting pressure evo- 
lution of the electronic state of CdV204 from orbital- 
ordered-insulator to orbital-ordered-metal and then to 
orbital-disordered-metal with increasing pressure in the 
tetragonal structure. 

As opposed to the electronic state, the magnetic state 
of the compound is found to be robust with increasing 
pressure as the MM of the V ion remains almost the 
same even at 105.5 GPa. In order to know the effect of 
pressure on the magnetic transition temperature we have 
plotted the Heisenberg exchange parameter J of neigh- 
bouring V ions as a function of pressure in Fig. 7. The 
pressure dependence values of J have been estimated by 
mapping of energies of FM and AFM solutions to the 
Heisenberg Hamiltonian and assuming that only neigh- 
bouring V ions have strong magnetic interactions. Keep- 
ing the localized nature of V 3d orbitals in mind it may be 
considered as a fairly good approximation. The J shows 
almost linear pressure dependence and it increases with 
increasing pressure indicating the enhanced exchange in- 
teraction strength due to increased overlap of V 3d or- 
bitals. If we scale the J with experimental Tjv then J at 
zero pressure is equivalent to Tv at zero pressure. Us- 
ing this scaling one can plot the pressure dependence of 
Tjv (P) /Tjv (0). For the direct comparison with the exper- 
imental data shown in Fig. 3 of Ref. 16, we have shown 
the interpolated value of Tjv(P)/Tjv(0) up to 10 Kbar in 
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the inset of Fig. 7. Our result shows about 2.5% increase 
in T N per GPa (i.e. 10 Kbar) of pressure which is closer 
to the experimental observation.— Keeping the complex 
nature of magnetic transition in mind where spin fluc- 
tuations due to geometrical frustrated can not be fully 
ignored, the applicability of such a simple mean-field ap- 
proach in understanding the pressure dependent shift in 
the Tjv is quite interesting and demands further investi- 
gation to such similar systems. 

IV. CONCLUSIONS 

The effect of pressure on the electronic and magnetic 
properties of tetragonal CdV2 04 have been studied by 
using GGA+C7 formulation of the density functional the- 
ory. The present study shows the importance of on-site 



Coulomb interaction among the V 3d electrons in deriv- 
ing the insulating ground state of the compound. The 
spin-orbit coupling is weak and responsible for titling 
the total magnetic moment (^1.3 fis) of the V ion by 
~9.5 degree with the z-axis. The d xz and d yz orbitals 
are found to be mainly occupied at neighboring V sites. 
The electronic structure of the compound shows interest- 
ing evolution from orbital-ordered-insulator to orbital- 
ordered-metal to orbital-disordered-metal with increas- 
ing pressure. The pressure induced metal-insulator tran- 
sition is essentially arising due to broadening of lower 
and upper Hubbard bands. The present work clearly 
demonstrate the applicability of the mean-field theory in 
understanding the pressure dependent behaviour of the 
magnetic transition temperature in similar geometrically 
frustrated systems. 
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V. FIGURE CAPTIONS: 

FIG. 1. (Color online) Atomic arrangement in the 
primitive unit cell. The tetrahedron contains four V 
atoms denoted by VI and V2. VI and V2 are orbitally 
inequivalent sites. 

FIG. 2. (Color online) Energy vs volume plots of FM 
and AFM solutions obtained from GGA+C/ ({7=4 eV) 
calculations. 

FIG. 3. (Color online) Partial density of states (PDOS) 
for (a) Cd 4d, (b) V 3d, and (c) O 2p symmetric states 
obtained from FM solution. Zero in the energy axis in- 
dicates the Fermi level. 

FIG. 4. (Color online) The total density of states 
(TDOS) per formula unit (fu) obtained from AFM so- 
lutions at 0.23, 16.6, 34.5, 45.3, 87.5, and 105.5 GPa are 
shown in figures (a), (b), (c), (d), (e), and (f), respec- 
tively. Zero in the energy axis indicates the Fermi level. 

FIG. 5. (Color online) Dispersion relations along high 
symmetric directions of the Brillouin zone for AFM so- 
lutions obtained for (a) 0.23 GPa and (b) 34.5 GPa. 

FIG. 6. (Color online) The electronic occupancies of V 
d xz and d yz orbitals along with orbial polarization P or b = 
{dxz-dy Z ) / {d xz +d yz ) for VI sites is shown as a function of 
increasing pressure. The insulating and metallic regions 
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are separated by vertical dashed line. obtained from GGA+U (U=4 eV) calculations. The in- 

FIG. 7. (Color online) The Heisenberg exchange pa- set shows the pressure dependence of T N (P)/T N (0) ob- 
rameter (J) for neighbouring V ions at different pressure tained from the scaling of the J. See the text for details. 
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